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ABSTRACT 

We present an analysis of the Evershed effect observed with a resolution of 0"2. Using the new Swedish 
1-m Solar Telescope and its Littrow spectrograph, we scan a significant part of a sunspot penumbra. Spectra of 
the non-magnetic line Fe I 7090.4 A allows us to measure Doppler shifts without magnetic contamination. The 
observed line profiles are asymmetric. The Doppler shift depends on the part of the line used for measuring, 
indicating that the velocity structure of penumbrae remains unresolved even with our angular resolution. The 
observed line profiles are properly reproduced if two components with velocities between zero and several 
km s"' co-exist in the resolution elements. Using Doppler shifts at fixed line depths, we find a local correlation 
between upflows and bright structures, and downflows and dark structures. This association is not specific of 
the outer penumbra but it also occurs in the inner penumbra. The existence of such correlation was originally 
reported by Beckers & Schroter (1969), and it is suggestive of energy transport by convection in penumbrae. 
Subject headings: convection - Sun: magnetic fields - sunspots 



1. INTRODUCTION 

The spectral lines observed in sunspot penumbrae are both 
asymmetric and shifted in wavelength. This phenomenon, 
known as the Evershed effect, is commonly attributed to the 
presence of large and highly inclined plasma flows. Although 
! the e ffect was discovered a hundred years ago (E vershed 
. 1 19091) the true nature of these flows remains enigmatic, which 
hampers a full understanding of the penumbral physics. The 
observational properties of penumbrae and the Evershed ef- 
fect a re well covered in a number of recent rev iew papers 
(e.g., Thomas & WeissI Il992l llOOl ISolankil 1200 3), and we 
refer to them for an overview. This introduction is focused 
on a particular aspect of the effect directly connected to our 
observational work. 

Penumbrae are almost as bright as the quiet Sun. Since 
the transfer of energy by radiation is inefficient, either the 
penumbrae are transparent and therefore very shallow, or the 
energy radiated away by penumbrae is transported from be- 
low by some sort of convective process^ Shallow penum- 
brae require nearly horizontal magnetic fields, which are 
, not observed and which leaves us only the second option 
(fSolanki & Schmidt 1993). Therefore, the penumbrae of 
sunspots, with their filaments and flows, are likely the re- 
sult of convection taking place in highly inclined strong mag- 
netic fields. However, there is no agreement on how the con- 
vection operates, i.e., on how the hot plasma rises, releases 
energy, and returns to the sub-photospheric layers, and how 
these motions occur in a strong magnetic field constraining 
the free movement of the plasma. Several possibilities have 
been put forward in the literature, and we will mention a 
number of them as example. The penumbrae may contain 
convective rolls, where the magnetic field lines are rigidly 
transported by the plasm a motions, and w hich require very 
inclined magnetic fields (iDanielsonI 1 1 96 Ih . Depending on 
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^ Here and throughout the text we adopt the loose definition of convection 
by Nordlund (2000); "convection is the transport of energy by hot fluid mov- 
ing up wards and cold fluid moving downwards". ^SDruil 1987i : ISchmidt et al] 



the magnetic field inclination, con vection sets in a s convec- 
tive rolls or traveling waves (Hur lburtetani2000t) . Siphon 
flows along magnetic field lines have been proposed mainly 
for explaining the hori zontal component of the Evershed 
flow jM eyer & Schmidt [19681 iThomas & MontesinosI [19931: 
IWeiss et al..2004 ). However, such flows transport hot plasmas 
from the sub-photosphere, which become optically thin and 
radiate away internal energy. Individual magnetic fluxtubes 
that raise, cool off, and then fall represent another possibil- 
ity. This process is c alled interchange convection by movin g 
magnetic fluxtubes dSchmidtl ll99ir iJahn & Schmidt* fT994l) 
and it has been modeled in detail by iSchlichenmaie r et al] 
(ll998af R). A significant part of the convective energy trans- 
fer is not provided by the interchange motions, but by flows 
along field lines indirectly induc ed during the interchange 

). In the vein of the thin 



jSchhchenmaier & Solaiikil 12003 
penumbra models dSchmidt et al. 



1J86), but improved to get 



rid of the problems. Spr uit & Scharmen (12006) propose the 
intrusions of non-magnetic convective cells. 

None of the modes of convection mentioned above satisfy 
all the observational constraints and, simultaneously, none of 
the observations seem to be so deciding as to rule out par- 
ticular models once and for all. The main problem probably 
lies in the insufficient spatial resolution of our observational 
description of the phenomenon. Convection could be easily 
identified upon detecting a clear correlation between vertical 
velocity and brightness, similar to that characteristic of the 
solar granulation. However, the flows in the penumbra are 
known to change direction and speed within scales smaller 
than the photon mean-free-path, a property inferred directly 
from observations (e. g., from the existence of broad-band 
circular polarization, ISanchez Almeida & LitesI Il992t see 
also § O. The photospheric mean-free-path length scale, 
some 150 km, is comparable with the angular resolution of 
the best images, therefore, even our best measurements are 
expected to average several structures with different physical 
properties. The observations always provide ill-defined 
averages of the measured properties, being each resolution 
element a volume of the solar photosphere. The size in the 
plane perpendicular to the line-of-sight (LOS) is set by the 
angular resolution, whereas the size alongthe LOS is set 



by the radiative transfer smearing (e.g., ISanchez Almeidal 
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Il998l) . Unscrambling the averages to retrieve properties of the 
unresolved structure is both difficult and prone to misinterpre- 
tations. Researchers have approached the problem assuming 
the exi stence of unresolved structures during the data analysis 
(e.g., LB umba 1960; Grigorjev & Katz 1972; Golovko 1974; 
Sanchez Almeida & Lites 1992; Solanki & Montavon 
1993; Wiehr 119951; ISanchez Alm eida et al. 1996; 
Westendorp Pla za et all 120011; Eellot Rubio, ,2004i: 
Sanchez Almeid^ 20051) . Then the outcome of the anal- 
ysis depends on the specific assumptions, a difficulty that 
increases as the resolution worsens and the observational 
average becomes coarser. On the contrary, improving the 
angular resolution removes part of the ambiguities and 
facilitates the interpretation of the observations. In view of 
past experiences, it is difficult to specify an observational 
plan that secures sorting out the penumbral puzzle. However, 
any successful plan necessarily requires improving the 
angular resolution to a point where the interpretation of the 
observations is not longer ambiguous. 

The advent of the Swedish 1-m Solar Telescope (SST; 
IScharmer et alj ^2003a b) has opened up new possibilities 
since it routinely provides a n angular resolution s ignificantly 
better than its predecessors (IScharmer et al.ll20021) . We take 
advantage of the SST to carry out spectroscopic observa- 
tions of the Evershed effect with an exceptional spatial res- 
olution (~0."2; § |2|i. Our study is focused on the lo- 
cal correlation betw een brightness and Doppler shift origi- 
nally found by Beckers & Schrotei^ ("1969) and later observed 
by others ( Sanchez Almeida et al. 1993; Johannesso n 1 1 993l 
Schlichenmaier & S chmidtli 1999c iSchmidt & Schlichenmaierl 
2000h . Features brighter than the local mean are associated 
with blueshift, and vice versa. The same correlation exist both 
in the limb-side penumbra and the c enter- side penumbra, a 
fact invoked by iBeckers & Schrote^ (1 19691) to conclude that 
vertical motions rather than horizontal motions are responsi- 
ble for the correlation. Such correlation between vertical ve- 
locity and intensity is characteristic of the granulation. The 
fact that the same correlation is also present in penumbrae 
suggests a common origin for the two phenomena, namely, 
convection. Despite the qualitative agreement, the amplitude 
of the velocity fluctuations, of the order of 100 m s~', is insuf- 
ficient to transport the energy radiated away by penumbrae. 
Velocities similar to those of the non-magnetic granulation 
are required independently of the specific mode of co nvec- 
tion responsible for the transport (see, e.g..'Spruill ll987i) . We 
revisit the issue with the best resolution available at present 
to (a) confirm the existence of such correlation, and (b) see 
whether the improved resolution reveals vertical velocities 
large enough to account for the radiative losses of penum- 
brae. With this in mind, we seek independent support for var- 
ious radically different observational studies suggesting the 
presence of convective motions in penumbrae. Filament like 
struct ures in velocity maps are found to join bright and dark 
knots dSchlichenmaier et al.ll2004il2005h . A complex velocity 
pattern with small-scale upflows and downflows is able to re- 
produce the asy mmetries of the Stokes profiles characteristic 
of penumbrae (Sanc hez Almeidall2005h . The same scenario is 
qualitatively consistent with the prop er motions observed in a 
penumbra bv lMarquez et alj (l2006bl) . The proper motions are 
predominantly radial but with a sideways component which 
arranges the plasma forming long narrow filaments co-spatial 
with dark features. Mass conservation arguments indicate that 
the accumulation of plasma in the dark filaments must be bal- 
anced by downflows of the order of 200 m s"' . 



Our study of the Evershed effect begins by describing the 
observations and the data reduction (§|2]and Appendix lAl). 
The shape of the observed line bisectors is analyzed and dis- 
cussed in § [3] The observed spectra show the correlation 
between Doppler shift and brightness discussed above. Its 
variation with the azimuth in the penumbra allows us to sep- 
arate vertical and horizontal velocities (§ |4|i. Line asymme- 
tries are analyzed in §|5]in terms of two spatially unresolved 
components. The implications of our work for understand- 
ing the structure of penumbrae are discussed in §|6l where we 
also consider consistencies and inconsistencies with previous 
studies. Finally, Appendix |B] shows that the kind of weak 
correlation between vertical motions and brightness found in 
penumbra is also observed in the quiet Sun when the spatial 
resolution is insufficient to resolve individual granules. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out on October 16, 2004, us- 
ing the SST equipped with AO and a Littrow spectrograph 
(Schai-meretal. 2003a b; Kiselman et al. 2006, We ob- 
tained spectrograms of the non-magnetic line Fe 1 7090.4 A in 
the penumbra of the decaying circular sunspot NOAA 10682 
(heliocentric angle 6 ~ 35°; cos6'=0.82). The black rectangu- 
lar box in Figure[T]indicates the field-of-view (FOV). We em- 
ploy a non-magnetic line to avoid contaminating the inferred 
velocity fields with cross-talks coming from the penumbral 
magnetic field. The penumbra was scanned from the inner 
to the outer boundaries with steps of 0."12 perpendicular to 
the slit. The slit width was set to 0"1 1, which corresponds to 
29 mA in the focal plane of the spectrograph and sets the spec- 
tral resolution of the data. Given the pixel size of the camera 
and the effective focal length at the focal plane of the spec- 
trograph, the sampling interval of the spectra is 0"04 along 
the slit and 10.5 mA across the slit, which suffices for the 
spatial and spectral resolution of the data. We use frame se- 
lection triggered by the spectrograph camera, providing the 
four spectra of highest contrast every 20.8 s. The exposure 
time was set to 150 ms. The scan covers an effective FOV of 
10."5 (across the slit) x 34."4 (along the slit), and it required 
3 1 min to be completed. Together with the spectrograms, we 
gather simultaneous slit-jaw and Ca H images (an example of 
the former is shown in Fig. [1]). They were used to set up the 
polar coordinate system used for the analysis of velocities car- 
ried out in § m which requires identifying the direction of the 
solar disk center, as well as assigning a curvature radius and 
a curvature center to the penumbra. The data are corrected 
for dark current and flatfield using standard procedures. In 
addition, we correct the spectra for the Modulation Transfer 
Function (MTF) characteristic of an ideal 1-m telescope at the 
working wavelength. Since the penumbral structures are pre- 
dominantly elongated across the slit, the deconvolution pro- 
cess can be reduced to a ID problem using as MTF a radial 
cut across the 2D axi -symmetric MTF of a diffraction h mited 
telescope (see, e.g., ISanchez Almeida & Bonetlll998h . The 
deconvolution, including optimum filtering, is carried out in 
the Fo urier domain followi ng the procedure explained in de- 
tail by Sobotka et all (Il993h . 

From the best spectrogram per slit position, we construct 
the 2D continuum image of the FOV shown in Figure |2^. 
Each column corresponds to one slit position so that the full 
set provides an image. Residual image motion in the direction 

TRIPPEL spectrograph ; see |http://www.solarphysics.kva.se/| 
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Fig. 1 . — Slit-jaw of one of the exposures showing our field of view (the box). The axes are given in Mm from the lower left comer, and an aiTow points out 
the direction of the solar disk center. The scan was carried out from left to right with the slit in the vertical direction. (The slit was artificially removed from this 
slit-jaw image during the flatfielding procedure.) The image was obtained with a 7 A wide color filter centered in 7057 A. 
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Fig. 2. — (a) Map of continuum intensity, (b) Doppler shift of the bisector at 80% of the line depth. Running means have been subtracted in these two maps to 
bring out only local variations, (c) Difference between the Doppler shifts of the bisectors at 80% and 20% without running mean subtracted. Images (b) and (c) 
shai'e the same scale of grays which goes from -2 km s"' (black) to +2 km s"' (white), (d) Original map of the Doppler shift without the running mean subtracted 
and represented with a scale spanning from -4 km s"' (black) to +4 km s"' (white). According to our convention, redshifts con'espond to positive Doppler shifts. 
The centers of the symbols in (a) indicate the positions of the representative line profiles shown in Figure|4] 
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Fig. 3. — (a) Power spectrum for the intensity variations along tlie slit. 
Tlie dotted line shows the average among all our penumbral spectra before 
restoration. The solid line represents the same average after restoration. The 
horizontal dotted line corresponds to the level of power at the cutoff imposed 
by the telescope diameter and, therefore, it indicates the level of noise. There 
is power in the original unprocessed image even beyond 5 arcsec"' , corre- 
sponding to a scale of 0"2. The restoration includes a noise filter that damps 
down the spectrum starting at this threshold, (b) The solid line shows ob- 
servations of the intensity variation across a representative narrow structure 
in one of the spectra. The symbols con'espond to a least squares fit of this 
profile with a Gaussian plus a parabolic background, the latter represented by 
the dashed line. The Gaussian has Of 19 FWHM, indicating that features of 
this width can be resolved in our spectra. 

parallel to the slit has been removed from the image by cross- 
correlation techniques. We shift the different columns so that 
the cross-correlation between adjacent columns is maximum. 
The same small re-arrange of the columns inferred from the 
continuum intensity is applied to all the maps displayed and 
discussed in the paper. This manipulation was carried out 
mainly for aesthetic reasons, since the residual image mo- 
tion has very limited impact on the results of our analysis, 
which is based on correlations between parameters obtained 
in the same pixels. The angular resolution of our spectra turns 
out to be very close to 0"2, as inferred from the cutoff fre- 
quency of the power spectra of continuum fluctuations in the 
spatial direction of the spectrograms. Figure [3^ shows the 
power spectrum for the variations of intensity along the slit. 
We average the power spectrum of all the columns forming 
the continuum image in Figure |2^. Both the original spec- 
tra and the MTF reconstructed spectra show signals above the 
noise at 5 arcsec"', corresponding to a scale of 0."2. The same 
conclusion can be drawn from Gaussian fits to the narrowest 
continuum features in the spectra. Figure [Sji shows one of 
these features reproduced with a Gaussian having a FWHM 
of0"19. 

The work aims at analyzing absolute Doppler shifts, which 
requires a careful calibration of the wavelength scale. First, 
the relative wavelength scale along the spectrograph slit was 
set using spectrograms taken while the telescope pointing was 
moving across the solar disk. The spectral line must have 
the same wavelength in these flatfield images, an informa- 
tion used to bring to a common scale all the spectra taken 
in each position of the scan. As we argue below, the spec- 
trograph was stable enough during the scan to allow a direct 
comparison of the wavelengths taken at different positions, 
which provides a consistent relative wavelength scale for the 
whole scan. Therefore, one can set a common absolute scale 
by knowing the absolute wavelength of only one spectrum. 
We set such scale assigning zero velocity to the mean line shift 
obtained from the minimum of Fe I 7090 . 4 A in the darkest 
parts of the umbra. According to BeckersI (Il977l) . such abso- 
lute velocity is uncertain by less than 100 m s~'. Umbral os- 
cillations may have an influence on the velocity scale thus ob- 
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Fig. 4. — Representative line profiles showing large asymmetries and the 
corresponding bisectors. Top: profiles of bright features. Bottom: profiles of 
dark features. Left: profiles on the limb-side penumbra. Right: profiles on the 
center-side penumbra. The actual positions are indicated by the symbols in 
Figure|2^. Wavelengths are in km s"' in our absolute velocity scale, whereas 
the intensities are given in units of the mean quiet Sun continuum intensity. 

taine d, but their a mplitudes are again smaller than 100 m s~' 
(e.g.. lLiteslll992b . The effect of the spatial stray-light is an- 
alyzed in Appendix lAl and it turns out to be negligible. We 
can assure the stability of the spectrograph during the scan be- 
cause of the position the H2O telluric line at 7094.05 A. The 
wavelength of the minimum of this line in the absolute scale 
used for Fe I 7090.4 A is constant with a rms fluctuation of 
only 90 m s"' . Most of these fluctuations result from uncer- 
tainties in determining the minima of the weak telluric line 
(equivalent width of only 5 mA). In short, our absolute wave- 
length scale provides velocities with an accuracy of the order 
of lOOms-i. 

3. BISECTORS AND UNRESOLVED VELOCITIES 

Spectral lines formed in atmospheres of constant velocity 
are symmetric and shifted as a whole, therefore, one should 
get the same Doppler shift independently of the part of the 
line used for measuring. Since this is not the case of our 
penumbral spectra, we conclude that even with 0."2 angular 
resolution a significant part of the penumbral velocity struc- 
ture remains unresolved. We bring up this fact by computing 
the bisectors of the intensity profiles. The bisector of an inten- 
sity profile is the line joining the mid points of the profile at 
each intensity level. (The intensity level is usually quantified 
using a percentage, with 0% representing the line core and 
100% the continuum.) Figure |4] shows a few representative 
examples of observed line profiles and bisectors. The bisector 
must be a vertical straight line if the atmosphere has a constant 
velocity, however, the bisectors of Fe I 7090.4 A are curved 
and inclined. They show different velocities at different line 
depths, indicating that the velocity structure is finner than our 
resolution. The lack of enough resolution is not an exception 
attributable to a few pixels. Figure. |2}; shows a map of the 
velocity difference between the bisector at line wings (80%) 
and the bisector at line core (20%). The map of differences 
follows a pattern with the characteristic filamentary appear- 
ance of the penumbra. The typical values differ from zero, 
and values as large as one km s"' are not unusual. The disper- 
sion of velocities responsible for these bisectors must be large, 
a conclusion that we try to quantify with the two-Gaussian fits 
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presented in §|5] 

As it is pointed out in § [T] our resolution element is actu- 
ally a volume, with the size in the plane perpendicular to the 
LOS set by the angular resolution, and the size along the LOS 
set by radiative transfer smearing. Consequently, observing 
slanted bisectors do not allow us to know if the gradients of 
velocity occur along or across the LOS. The two of them prob- 
ably coexist (see §|6ll. 

4. RELATIONSHIP BETWEEN CONTINUUM 
INTENSITY AND VELOCITY 



In agr eement with iBeckers & Schroted (Il969h and oth- 
ers (e.g. Sanchez Almeida et al.' 'I 993I: iJohannessonI 1 1 9931; 
[Schmidt & Schlichenmaier 2000, see § [TJ, we find a local 
correlation between continuum intensity / and Doppler shift 
Ud- The purpose of the section is twofold. First, it describes 
the observed correlation. Second, horizontal and vertical mo- 
tions are separated to show that a significant part of the corre- 
lation is due to vertical motions. 

The local variations are obtained by subtracting from the 
original maps a running mean. Figure |2]3 has been computed 
from the original Figure |2}i using a running box 1 ."6 wide. 
The same box is used for the intensity in Figure |2^. We also 
tried with boxes half and twice this value to conclude that the 
actual width is not critical (see below). The correlation can be 
observed directly from inspection of the velocity and inten- 
sity images, e.g.. Figure |5] It shows how local blueshifts are 
preferentially associated with local bright features. The corre- 
lation is more clear in the limb-side penumbra (Fig.lSj"), but it 
is also present with the same sign in the center-side penumbra 
(Fig.|5^). The relationship can be quantified as, 



UD-{UD)^m(I-{I)) + k, 



(1) 



with the angle brackets denoting local running mean averages. 
In this equation and throughout the paper the intensity / is 
referred to the quiet Sun continuum intensity and therefore is 
a dimensionless parameter, whereas the sign of Ud is chosen 
to be positive for redshifts. 

Figure |6] shows the observed variation of the Doppler shift 
versus continuum intensity for two sections of our FOV, one 
in the limb-side penumbra (a) and the other in the center-side 
penumbra (b). We represent the mean Doppler shift con- 
sidering intensity bins of 0.01. The slope m varies within 
the sunspot, from the center-side penumbra to the limb-side 
penumbra. However, it maintains a negative sign implying 
that bright features are associated with blueshifts (Fig. |6]l. 
The fact that m keeps the sign implies that the correla- 
tion is produced by vertical motions, since the line-of-sight 
component of the horizontal radial velocities changes sign 
from the center-side penumbra to the limb-side penumbra 
((Beckers & Schroter 1969). 

This qualitative argument can be quantified adopting the 
following model for the relationship between velocities and 
intensities. 



U: = m,{I-{I)) + k„ 
Uh=mh {I- {I}) + kh, 



(2) 



where ;«,, k,, m/, and k/, are constants, and t/, and Ui, represent 
the vertical and horizontal velocities in the absolute reference 
system portrayed in Figure |7] {Ur>Q for upward velocities 
and Uh > for outward velocities.) If the motions are pre- 
dominately radial so that the azimuthal component can be ne- 
glected, 

UD = -(UhSmecos^ + U-^cose), (3) 



with the s ymbol 9 representing the hel iocentric angle of ob- 
servation (IMaltbylfT964t ISchrotedl 1 967h . This model leads to 
equation ([T]l with the slope m. 



m(0) = -(nif, sin 6 cos (f)+m^ cos 9) , 



(4) 



and with the running mean average Doppler shift, 

(Ud) =-{kh sin 9 cos (p+k^ cos 9). (5) 

The average Doppler shifts and slopes depend on the azimuth 
(p. In order to apply the model to the observed velocities, we 
construct a polar coordinate system that assigns azimuths and 
radial distances to each position of the FOV. The center of cur- 
vature of the observed penumbra was estimated visually by 
trial-and-error so that the grid of polar coordinates provides 
radii parallel to the observed penumbral filaments. The best 
compromise is shown in Figure |8] The direction of the solar 
disk center sets the origin of the azimuths and it was obtained 
by comparison of the slit-jaw images with quasi-simultaneous 
full disk MDI continuum images. Figure |9] shows m(0) ob- 
tained within the sectors of approximately constant azimuth 
in Figure |8] when the Doppler shifts in Figure ^ are used 
(bisector shift at line wings). Figures and|6]3 illustrate the 
kind of fit leading to m{(p). Using the model in equation dU, a 
linear fit m((/)) versus cos (p renders. 



m,~ (3.3 ±0.4) km s" 
mA~(-3.7±0.8)kms" 



(6) 



The error bars are 1 -sigma formal errors provided by the least- 
squares routine, and assume that the deviations of the data 
from the fitted line are due to random noise (see Fig.|9]). The 
sign mi, < indicates that horizontal velocities are enhanced 
in dark features, whereas m, > implies that upflows occur 
in bright features. The fact that ;«/, < and > is a firm 
result, as it is indicated by the error bars. We also checked 
that bisectors at different line depths can be used to compute 
velocities without modifying the resulting signs. The use of 
bisectors closer to the line core only reduces the amplitude of 
the velocities. We also repeated the calculations using only 
the inner half of the penumbra, i.e., the half closest to the 
umbra. The signs of m, and m/, are not modified. In addition, 
we carried out the whole analysis doubling and halving the 
width of the box used to compute the running mean. Again 
the signs of m, and m/, remain. Decreasing the width reduces 
the amplitude of the correlation whereas increasing the width 
has almost no effect on and m/, . 

With the same approach leading to and m/,, one can use 
the local average Doppler shifts (Ud) to infer the constants k^ 
and kh (see equation ||5l). It renders, 



A;, = (-0.48 ±0.05) km s" 
A:/, = (3.30±0.15)kms"' 



(7) 



where the error bars are also formal errors from the lin- 
ear least-squares fit. The finding of downward velocities in 
penumbrae (k^ < 0) is not new. Downward velocities have 
be en measured b y many o bservers starting with the work 
bv lRimmelell[T995l: see also IServajeanl [19611) . These down- 
ward motions tend to be localized in the external parts of 
the penumbra, but the actual distribution depends on the spe- 
cific measureme nt. Figure [TOh shows several among the 
pub lished values ([Rim mele 1995; Schlichenmaier & Schmidt 
[20001: iTritschler et al.1 12004t ISanchez Cuberes et alJ 120051) . 
which are presented as vertical velocity versus distance to the 
sunspot center. In order to compare our result with them, we 
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Fig. 5. — Continuum images of a portion of the center-side penumbra (a) and tlie limb-side penumbra (b) overlaid with contours of the velocity field inferred 
from bisector shifts at line wings (Fig. |2). Both the intensities and the velocities represent local variations; the large scale variation has been removed by 
subtracting a smoothed version of the original maps. Blue and red contours correspond to blueshifts and redshifts, respectively (±250 m s"'). The velocity 
contours show a clear filamentary structure. Note that blueshifts tend to be associated with bright features, and vice-versa. It is not a one-to-one correlation, 
though. Such association is more pronounced in the limb-side penumbra, but it is also present in the center-side penumbra, indicating that it must be produced 
by vertical rather than horizontal velocities. 



consider that and ki, in equation (|2]l depend on the radial 
distance to the sunspot center p. Our penumbra is divided 
in bands of constant radial distance (Fig. |8]l. Then equa- 
tion (|5]l was fitted to the azimuthal variation of the running 
mean Doppler shift in each one of these bands, which ren- 
ders k-{p) and ki,(p). The results have been represented in 
Figure [TOj the solid lines with symbols, with the symbols in- 
dicating the center of the band used to carry out the fits. We 
find mild upflows in the inner penumbra and downflows in 
the outer penumbra. These features are common to almost all 
the other measurements portrayed in Figure [TOk . The mean 
horizontal velocities reach values of up to 3.7 kms"' not far 
from the external penumbral border, in agreement with values 
reported in the literature (Fig.fTOb). Our curves remain within 
the scatter of the published values, a fact pointed out to argue 
that no obvious error burdens the validity of our velocity scale. 
The quantitative differences between the curves represented in 
Figure [To] can be attributed to many factors biasing this type 
of observation, e.g., variations from sunspot to sunspot, differ- 
ent systematic errors in the absolute velocity scales, the use of 
different spectral lines and procedures to infer velocities, the 
different ways of estimating sunspot radii, etc. Although the 
vertical velocities that we find are somehow larger than previ- 
ou s estimate s (Fig [TOk ) , the values are close to t hose r eported 
by lRimmeld (Il995h and lSanchez Cuberes et aP (l2005h . which 



correspond to sunspots observed at the very solar disk cen- 
ter, so that the Doppler shift directly yields vertical velocities 
without geomet rical transformat ions. In the particular case 
of Sanchez Cub eres et al.l (l2005h . they carry out a ID inver- 
sion of the spectra to interpret the observed line asymmetries. 
Consequently, they obtain a velocity for each optical depth 
along the LOS. We represent the velocity at the bottom of the 
atmosphere (continuum optical depth equals 0.3), where the 
line-wings are formed. 

We would like to stress that interpreting the average vertical 
velocities shown in Figure [TOh is not trivial. They have been 
inferred as systematic Doppler shifts of the observed spectral 
lines (the part of the Dopple r shift indep endent of cos (/>; see 
equation 0)- According to ISanchez Al meida (2005), these 
net Doppler shifts may not necessarily represent net vertical 
plasma motions. If the penumbrae have plasma moving up 
and down, and if the thermodynamic conditions of these up- 
flows and downflows are not identical, then there will be a net 
displacement of the average spectrum even though there is no 
net plasma motion. This lack of cancellation of the net sig- 
nal is responsible for the well known convecti ve blueshift of 
the sp ectral lines observed in the quiet Sun (e.g. lDravins et al.l 
[Ml . 

5. TWO-GAUSSIAN FIT 
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Fig. 6. — Doppler velocity versus continuum intensity in sectors at tlie 
limb-side penumbra (a) and the center-side penumbra (b). The symbols and 
the solid lines coirespond to observations and linear fits, respectively. The 
error bars represent standard deviations of the mean Doppler shift considering 
all the points with a given intensity /—(/). Velocities are measured from 
bisectors at line wings (Fig.|2j)). 
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Fig. 8. — Image of the FOV with the grid used to assign to the spectra polar 
coordinates centered on the sunspot umbra. 
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Fig. 7. — Schematic with a sunspot and the reference system used to define 
angles and velocities. The line-of-sight (LOS) corresponding to the heliocen- 
tric angle of observation 9 is shown as a dashed line. The velocity U of a 
point on the penumbra with an azimuth cf> only has a vertical component U- 
and a radial horizontal component Ui, . The sunspot axis coincides with the 
solar vertical direction, and the line joining the sunspot axis and the solar disk 
center sets the origin of </>. 



Fig. 9. — Slope of the relationship between intensity and Doppler shift 
m(<f>) versus the cosine of the azimuthal angle within the sunspot <!>. The error 
bars are formal error bars obtained from the fit used to derive m((l>). 



We also studied whether the observed Hne asymmetries can 
be reproduced by means of two spatially unresolved velocity 
components. Rumba ( 1960) puts forward the idea, and it has 
been pursued by many oth ers ever since (e.g., Maltbv 1964t 
IStellmacher & Wiehdll97UlWiehr 1995; .Ichimoto 1988) . We 
fit the line profiles 5(A) with two Gaussian functions plus a 
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Fig. 10. — (a) Variation of k- with the distance to the sunspot center (the 
solid line with symbols). Note the tendency for k- to be negative. It is not very 
different from other Doppler shifts reported in the literature and also included 
in the plot. The symbo ls in the inset stand for Schlichenmaier & SchmidJ 
EOGO) (S&S GO), Rimmele 1 1995) (R 95), iTnFschler et al.l 12004^ (T et al 
04), and Sanchez Cuberes et al. (2005) (SC et al 05). (b) Variation oik,, with 
the distance to the sunspot center All radial distances are refeired to the 
sunspot radius Rs (~14 Mm). 

continuum, 

5(A)=5,(A)-^!'iexp{-i[(A-Ai)/ai]2} (8) 
-fl'2exp{-i[(A-A2)/a2]2}, 
with the continuum given by, 

5t-(A) = ao + A ai + A^ 02, (9) 

and symbols A,-,cr,- and di representing, respectively, the cen- 
tral wavelength, the width, and the depth of the two Gaussians 
(/ = 1,2). The fits are in general very good, like the examples 
shown in Figure [TT] Note how the symbols used to represent 
the observations cannot be distinguished from the dotted line 
representing the fits, meaning that the two Gaussians are able 
to reproduce the asymmetries of the observed line profiles. 

The two components have very different Doppler shifts 
(i.e., different A,). The components in each pixel have been 
classified either as shifted or as unshifted depending on which 
one has the largest absolute Doppler shift. As it is shown 
in Figure [121 the unshifted component has almost no veloc- 
ity throughout the penumbra, whereas the shifted component 
follows the pattern characteristic of the Evershed effect, with 
blueshifts in the center-side penumbra and redshifts in the 
limb-side penumbra. The velocities of the shifted component 
are larger than those inferred from the bisectors. 

The two components contribute similarly to the line pro- 
files, except in the inner penumbra, where the unshifted com- 
ponent dominates. An example of two-Gaussian fit in the in- 
nermost penumbra is shown in Figure [TT] bottom. Note the 
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Fig. 1 1 . — Examples of two-Gaussian fits. Top: limb-side profile. Middle: 
center-side profile. Bottom: umbra-penumbra border profile. The symbols 
represent the observations whereas the dotted line shows the fits. The solid 
lines coiTespond to the two Gaussians used in the fit. The continuum has been 
added to the individual Gaussians to facilitate comparison. 

broad and shallow shifted component, and how the unshifted 
component dominates the line profile. Figure [13] illustrates 
the same result in a more systematic way. We use the equiv- 
alent width to parameterize the strength of the component. 
The equivalent width W is defined as the area above the line 
profile normalized to the continuum intensity. In the case of 
Gaussian profiles, 

Wi^V2^diai/S,(X,), (10) 

where we have neglected the variation of the continuum 
within the line profile so that the continuum intensity SdX) ~ 
SciXc), with Af the line core wavelength. Figure [13] shows 
the excess of equivalent width of the two components as a 
function of the distance to the sunspot center Following the 
approach of the previous sections, we remove a running mean 
local average to show local variations. We subtract {W)/2 
with (W) = (Wi) + {W2)- This local mean is the same for the 
two components to ensure that the differences do not come 
from the two components having different local means. Fig- 
ure [13] shows how the unshifted component tends to have an 
equivalent width significantly larger than the shifted compo- 
nent for p < 0.65 Rs, whereas the contribution of the two com- 
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Fig. 12. — Mean Doppler shift versus cosine of the azimuthal angle for the 
two Gaussian components that reproduce the line asymmetries, (a) Shifted 
component. It shows the pattern characteristic of the Evershed flow, with 
blueshifts in the center-side penumbra (cosiji ~ 1) and redshifts in the limb- 
side penumbra (cos(/> ~ -1). (b) Unshifted component, showing very small 
Doppler shifts. Each point represents the average Doppler shift over a small 
range of azimuths and radial distances (the sectors represented in Fig. [S). 
Error bars are not included to avoid overcrowding, but they are smaller than 
the scatter among individual points. 



ponents balances from this radius on. (The symbol Rs stands 
for radius of the sunspot.) Actually, the shifted component 
shows a small excess at 0.75, and then the unshifted compo- 
nent takes over again. 

The unshifted component seems to be brighter since it dom- 
inates in the bright pixels. Figure[T4]shows how the equivalent 
widths of the two components vary with the local intensity. 
The equivalent width of the unshifted component is larger in 
the bright pixels (/ > (/)), whereas the contribution of the two 
components balance in the locally dark features (/ < (/)). 

As we point out above, the Doppler shifts of the shifted 
component resembles but exceeds the Doppler shifts derived 
from the bisectors. The corresponding mean vertical and 
horizontal velocities can be inferred following the same ap- 
proach employed for the velocities based on bisectors. Fig- 
ure [15] shows k- and ki, obtained from a least squares fit of 
the Doppler shifts of the two components using equation (|5]l. 
Points grouped by radial distance are used to obtain the radial 
variation. (The same approach used to derive Fig. [TO]) The 
shifted component shows systematic horizontal velocities of 
up to 5 km s"', and it presents moderate upflows in the inner 
penumbra {k- < 0.5 km s"' ) and larger downflows in the outer 
penumbra (k~ ~ 1.5 km s"'). Simultaneously, the unshifted 



Fig. 13. — Radial variation of the equivalent widths of the shifted com- 
ponent (the dashed line) and the unshifted component (the dotted line). The 
same local mean equivalent width has been subtracted from the two compo- 
nents. 
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Fig. 14. — Equivalent width of the shifted component (the dashed line) and 
the unshifted component (the dotted line). The unshifted component dom- 
inates in locally bright structures, suggesting an association between bright 
features and the unshifted component. The error bars represent the standard 
deviation of the averages represented in the figure. 



component shows velocities never exceeding 1 km s"'. Be- 
cause of the quality of the fits, one can conclude that repro- 
ducing the observed bisectors requires plasmas with velocity 
amplitudes between and 5 km s"' coexisting in our 0."2 res- 
olution elements. In the case of the shifted component the 
upflows in the inner penumbra exceed our observational un- 
certainties (^100 m s"'; see §|2]i. However, we are reluctant 
to interpret these shifts as systematic upflows present in the in- 
ner penumbra. If the properties of the shifted component are 
taken literally, they indicate the existence of spatially unre- 
solved upflows and downflows of the order of several km s"', 
a feature that is more meaningful that the small systematic up- 
flow. The upflow may be real, but it may also be a false resid- 
ual left if the Doppler signals of the large upflows and down- 
flows are not strictly proportional to the flow of mass. The 
existence of such large upflows and downflows follows from 
the analysis of Figure [16] which shows the variation with the 
radial distance of the widths of the two Gaussian components 
(cTi and (72 in equation [SI). The width of the shifted compo- 
nent reaches a value of up to 5 km s"' in the inner penumbra, 
and this value is almost independent of the azimuthal angle 
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Fig. 15. — (a) Vertical velocities k- and (b) horizontal velocities kj, for the 
two components resulting from the two-Gaussian fits of the asymmetric line 
profiles. We also include the velocities inferred f rom the bisector analysis for 
reference (the solid line with symbols, as in Fig. [To]) The unshifted compo- 
nent presents very small velocities. Systematic horizontal velocities of up to 
5 km s"' characterize the shifted component (the dashed line). This compo- 
nent also shows upflows in the inner penumbra (k- > for p/Rs < 0.7) and 
downflows in the external penumbra. 

within the sunspot. This width is too large to be due to thermal 
motions (~ 1 km s"' for penumbral temperatures of 5500 K), 
and it is also much larger than the typical widths of the pro- 
file of Fe 1 7090.4 A (-2 km s"' in the quiet Sun). Then if 
the inferred widths are due to spatially unresolved motions, it 
means that redward and blueward velocities of several km per 
second are associated with the shifted component in the inner 
penumbra. Since the widths do not depend on the azimuthal 
angle within the sunspot, such large velocity dispersion is also 
representative of the vertical velocities and, consequently, up- 
flows and downflows of several km per second are associ- 
ated with the shifted component. According to Figure [161 
the width of the shifted component suffers a strong decrease 
from the inner to the outer penumbra. Observations have con- 
sistently shown that the widths of the sp ectral lines increase 
from the inner to the outer penumbra ( e.g., |johannessonll993l : 
lRimmelel[T99l iTritschler 6711112004 . which seems to be at 
variance with our result. However, there is no obvious incon- 
sistency keeping in mind that the shifted component does not 
contribute to the line profile in the inner penumbra (Fig.[T3]l. 
The profile is dominated by the unshifted component, which 
is very narrow and whose width increases with the radial dis- 
tance (Fig.[T6]l. Consequently, the line width of the full profile 
is expected to go from the values of the unshifted component 
in the inner penumbra, to some kind of average between the 
widths of the two components in the outer penumbra, thus 
increasing from the inner to the outer penumbra (see Fig. [TSI l. 
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Fig. 16. — Variation with the radial distance of the widths of the two Gaus- 
sian components. They are given in km s"' versus radial distance normalized 
the Sunspot radius R,. The dotted line corresponds to the unshifted compo- 
nent whereas the dashed line represents the shifted component. (A Gaussian 
has a FWHM some 2.35 times the parameter cr,- represented in the figure.) 

A final caveat is in order Like the use of bisectors, the two- 
Gaussian fit is only an heuristic representation of the observed 
line asymmetries. The two components resulting from our fits 
do not necessarily represent two physical components exist- 
ing in the penumbra, and a literal interpretation of the prop- 
erties of the components may be misleading. The Evershed 
effect is not only the result of combining two separate com- 
ponents. This description is incomplete since two unresolved 
but symmetric components cannot account for the broad-band 
circular polarization observed in sunspots (see §|6]l. 

6. DISCUSSION AND CONCLUSIONS 

We find local upflows associated with bright structures and 
local downflows associated with dark structures in 0."2 angu- 
lar resolution spectra (§ |4]i. The correlation keeps the same 
sign in the limb-side penumbra and the center-side penumbra, 
which indicates that it is due to vertical motions. The correla- 
tion remains even when only the inner half of the penumbra is 
considered. The existence of such corr elation was originally 
reported by iBeckers & Schroteil (1 19691) . and it is suggestive 
of convective motions in penumbrae. It has the sign required 
to transport energy by convection, although the magnitude of 
the inferred velocities is not yet enough to account for the ob- 
served radiative flux (see the discussion below). This deficit 
could be caused by the still insufficient spatial resolution of 
the present data. Appendix |B] shows how the correlation in 
penumbrae is qualitatively similar to the intensity velocity re- 
lationship found in the quiet Sun when the angular resolution 
is not enough to resolve individual convective cells (granules). 
In the case of the quiet Sun the correlation is produced by 
convective motions. The fact that the residual correlation is 
similar to that observed in penumbrae does not prove the con- 
vective nature of the penumbral pattern. However, the com- 
parison is very suggestive. The relationship between verti- 
cal velocity and intensity is not the only observable hinting at 
convection in penumbrae. For example, the proper motions 
arrange the plasma forming long narrow filaments co-spatial 
with dark features, in a behavior char acteristic of t he convec- 
tive proper motions of the quiet Sun (iMarquez et al. 2006b). 
Another hint comes from the observed Stokes asymmetries, 
which can be reproduced by small magnetic loops carrying 
mass flows that emerge and return to the sub-photosphere 
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within the penumbra dSanchez Almeidall2005h . If the loops 
connect a hot footpoint with a cold footpoint, then the flow 
along magnetic fie ld lines produce s a vertical transport of 
energy (e.g., Schlic henmaier & So lanki 2003). Such loops 
may be small scale counterparts to the filament like struc- 
tures in velocity map s found to join bright and dark knots by 
ISchhchenmaier et alj (l2005h . 

We would like to emphasize that any of the modes of con- 
vection described in § [T] is in principle consistent with the 
observed correlation between vertical velocity and intensity. 
The predictions of the various models are simply not realistic 
enough to allow a quantitative comparison with our observa- 
tions, which should be regarded as a constraint to be fulfilled 
by any further realistic modeling. 

According to equations and (|6]), the relationship be- 
tween the continuum intensity fluctuation A/rms and the fluc- 
tuations of vertical velocity At/- are, 

Af/, ~3.3kms"'A/,™. (H) 

We observe A/rms to be of the order of 0.09, rendering 
a AUr ~ 0.3 km s"' insufficient to account for the radia- 
tive losses of penumbrae. Vertical velociti e s of t he order 
of 1 km s"' are required; see, e.g., IS pruitI (|1987|) . How- 
ever the fluctuations of continuum intensity are underesti- 
mated due to the insufficient resolution. If the true fluc- 
tuations are as large as 0.2 rather than the observed value, 
equation (fTTT i yields AU~ ~ 0.7 km s~', which could eas- 
ily account for the radiative losses of penumbrae. Such 
large continuum contrast is characteristic of numerical sim- 
ulations of non-magnetic convection in the quiet Sun (e.g., 
Utein & Nordlund 1998 obtain A/ms ^ 0.16). The above esti- 
mate implicitly assumes the solar surface to be fully covered 
by upflows and downflows. This simplifying hypot hesis is 
not ne e ded. Following arguments similar to those by ISpruitI 
(ll987h . lSch lichenmaier & Solanki ("2003) estimate the energy 
transpor ted b y the moving fluxtube model of interchange con- 
vection (ISchlichenmaier et al.lll998al) . The model fluxtubes 
have very fast hot stationary upflows which are shown to be 
capable of explaining the penumbral radiative losses even if 
the fluxtubes fill only a small fraction of the surface (a cov- 
erage of so me 10% suffices for upflow s of 4 km s~'). The 
estimate bv ISchlichenmaier & SolankH (12003) is also valid in 
a general case of concentrated upflows and downflows. They 
can be consistent with the radiative losses of penumbra and 
with our observation, provided that the pattern of concen- 
trated upflows and downflows, once smeared to our spatial 
resolution, gives rise to the intensity vertical velocity correla- 
tion described in § |4] 

The velocities discussed above are computed using the 
shift of the bisector in the wings of the non-magnetic line 
Fe I 7090.4 A . The use of bisectors closer to the line core 
does not modify the scenario described above, although the 
velocities are smaller. However, the fact that the bisector 
shift depends on the intensity level within the line is con- 
sequential. It proves in a direct and robust manner that the 
penumbral velocity structure remains unresolved in our spec- 
tra. This result is not unexpected, though. Images of penum- 
bra show structures at least a factor two thinner than our 
resolu tion (IScharmer et al.ll2002t iRouppe van der Voort et al.l 
|2004|) . and it is conceivable that part of such structuring is 
shared by the velocity pattern. Even more direct evidence for 
unresolved velocities comes from the presence of broad-band 
circular pola rization (BBCP) in penumbrae, discovered thirty 
years ago bv lllling et al.l (Il974allbl) . The BBCP demands large 



gradients of velocity along the line-of-sight within the range 
of heights where a typical photospheric fine is formed (e.g . , 
ISanchez Almeida & Liteslll992t ISolanki & Montavonlll993l) . 
This range is not larger than 150 km and therefore unless the 
spatial resolution is significantly better than this value, spec- 
tral lines with slanted bisectors are to be expected. 

Even our best spatial resolution does not allow us to re- 
solve the Evershed velocity pattern, which rises an obvious 
concern. All analyses of penumbral spectra assuming a sin- 
gle resolved component may lead to biased velocity estimates. 
The existence of a single component is an implicit assump- 
tion of all those measurements that assign a single velocity 
to each resolution element. Often this simplifying assump- 
tion allows to observe with a spatial re solution otherwise im- 
possible to reach (e.g. lLanghans et £11 12005). However, these 
measurements are liable to a significant bias since reproduc- 
ing the observed line shapes requires plasmas with veloci- 
ties between and 5 km s"' co-existing in the resolution ele- 
ments. This range of velocities is deduced from the two com- 
ponent fits carried out in § |5] but it is also very consistent 
with the dispersion of velocities inferred fro m the asymme- 
tries of the Stokes profiles (e.g., Westendorp P laza et al.l2001l; 
iMathew et £0120031: ISanchez Almeida 2005.) . If the proper- 
ties of the observed spectra are measured assuming a single 
component per resolution element, the measurement only pro- 
vides an ill-defined average of the unresolved velocities. Dis- 
crepancies between the measured velocities and the true mass 
weighted average velocities are to be expected, with differ- 
ences that can be a fraction of the dispersion of velocities in 
the resolution element. Since the dispersion of velocities in 
the penumbra is of the order of a few km s"', differences as 
large as one km s"' are conceivable. Obviously, this caveat 
also affects the velocities measured when working out the 
intensity velocity relationship in § |4] Note, however, that 
these are differential measurements unaffected by a global 
bias. Only differences of the bias between bright and dark 
features are of concern, and they are expected to be of second 
order. 

In addition to the correlation between intensity and verti- 
cal velocity pointed out above, we find two other results in 
good agreement with previous observations. First, the hor- 
izontal velocity increases in dark penumbral features. This 
conclusion follows from the negative slope of the correlation 
between horizontal motions and intensities (m/, < in equa- 
tion 111). Second, we find moderate but systematic down- 
flows in the penumbra^. The presence of such downflows 
over a significant part of the penumbra is a result consistently 
shown by most recent measurements (see Fig. [Toh and the 
references in there). We want to emphasize these two non- 
trivial coincidences with previous works before discussing ap- 
parent disagreements with two other recent observational pa- 
pers. Such disagreement is particularly discomforting since 
the works also employ high spatial resolution SST data. Be- 
cause of the potentia l importance, w e discu ss the disagree- 
ment in some detail. iLanghans et al.l ( 120051) do not find the 
systematic downflows in penumbrae discussed above. Their 
measurements are optimized for magnetic studies, therefore, 
they employ a magnetic line whose changes of shape due to 
magnetic field variations may induce spurious velocity sig- 
nals. They use filtergrams with moderate spectral resolution 

As we argue in the previous paragraph, this may not represent a net ver- 
tical motion of the plasma but a bias produced by the imperfect cancellation 
of unresolved upflows and downflows. 
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(72 mA FWHM), thus averaging information coming from 
different line depths. Only two wavelengths within the line 
profile are sampled, and these are not far from the line core 
(±50 mA). These parameters suggest a weighting of their 
measurements towards the line core, where the systematic 
downflows that we find are milder (the bisector at the 20% in- 
tensity level has downflows only 1/3 of the values in Fig.fTOb). 
All these factors combined may eventually explain the lack 
of systematic downflows. Bellot Rubio et al. (2005) find dark 
cores in penumbral filaments having upflows instead of the 
downflows to be expected if a dark struct ure follows the in- 
tensity vertical velocity relationship in §|4] iBellot Rubio et al] 
( l2005h employ the same instrumentation used in our work. 
Even the same spectral line is used in some of the estimates. 
Consequently, the issue of spectral resolution argued above 
is of no relevance for sorting out discrepancies. The differ- 
ences must be pin down to the technique of analysis. One 
can identify two main differences with respect to our case. 
First, they use as absolute velocity reference the wavelength 
of the core of the spectral l ine in the dark cores. According 
to lBellot Rubio et al] ( l2005h . this reference provides an abso- 
lute reference within a "few hundred m s"' ". However, in the 
inner penumbra, where the dark cores reside, a shift of the ab- 
solute velocity scale by 200 m s"' suffices to turn upflows into 
downflows and vice-versa (see Fig.fTOh for p/Rs <0.7). Sec- 
ond, they study a limited number of specific dark cores. The 



correlation that we find is not one-to-one and those cases may 
correspond to deviations from the mean law. In this sense, the 
dark features giving rise to the correlation between intensity 
and velocity may not correspond to dark cores but a different 
penumbral dark structure. In short, the reasons for the ap- 
parent discrepancies between our results and these two SST 
based works are so far unclear and must be investigated with 
new observations. They will also allows to correct an obvious 
limitation of our work. It is based on one half of the penum- 
bra of a sunspot observed in a single position on the disk. 
The used of a data set with such poor statistics is justified by 
the exceptional angular resolution of the spectra, and by the 
fact that Evershed flow seems to be a property common to all 
penumbra. However, the results require confirmation. 
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APPENDIX 

EFFECT OF STRAY-LIGHT ON THE ABSOLUTE WAVELENGTH SCALE 

The spatial stray-light contaminates the average umbral spectrum used to set the absolute wavelength scale. The observed 
umbral spectrum S{\) is not the true umbral spectrum 5„(A) but it contains a fraction a of stray-light spectrum Ss{X), 



5(A) = a5,(A) + (l-a)5„(A). 



(Al) 



The contamination artificially shifts 5(A) with respect to Su{X) if ^^(A) is shifted with respect to 5„(A) by an amount As ^ 0. For 
the sake of simplicity, we assume the stray-light spectrum to have a shape similar to the umbral spectrum, therefore. 



5,(A)~C-'5„(A-A,), 



(A2) 



with the constant C standing for the ratio between the continuum intensities of the umbra and the stray-light. To first order in A 



5,(A) ~ C"' 

and inserting the previous expression into equation ( lAll i. 



5„(A)-A,- 



5(A)~(l-a + a/C)5„(A-Aa), 

An. = Ac 



1+ C 



a 



(A3) 

(A4) 
(A5) 



where Aq. gives to the artificial shift introduced by the stray-light. In order to estimate Aq, one needs to know the fraction of 
stray-light a, the shift of the stray-light contamination A,, and the continuum contrast of the stray-light C. We consider 



a < 0.05, 



(A6) 



to be a realistic upper limit to the stray-light contamination. The brightness of the Venus disk in SST images taken during the 
June 2004 Venus transit^ provides the level of SST stray-light. At the wavelength of our observation, the Venus disk shows an 
intensity smaller than 0.05 times the quiet Sun intensity as soon as it is measured 1" inside the Venus limb. The umbral stray-light 
signals come from the penumbra or beyond. The umbral region used for wavelength calibration is separated from the penumbral 
border by more than 1 ", which justifies the upper limit in equation ( IA6I ). The shift of the stray-light profile is unknown but, if the 
stray-light signals come from the neighboring penumbra, it can be estimated from the line core shift of the observed penumbral 
profiles. The shift of the mean penumbral profile turns out to be of only 30 m s"' . The shift considering only the inner penumbra 
is even smaller These figures are effected by the uncertainty that we are trying to estimate and, therefore, they must be taken 
with some caution. However, the mean penumbral profile making up the stray-light contamination is expected to present a shift 
much smaller than the large Evershed shifts discussed in the paper First, we are using the line core for wavelength calibration. 



* http://www.solarphysics.kva.se/. 
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where the Evershed shifts are largely reduced. Second, the Evershed shifts of the center-side penumbra and limb-side penumbra 
have opposite signs so that they tend to cancel in the averages. Only the (small) vertical component of the Evershed effect is left. 
We use as an upper limit for A., the convective blueshift of Fe I 7090.4 A in the quiet Sun, 

|A,| < 300 ms-'; (A7) 

see lDravins et al.l(ll981l Fig. 6b). By considering this value, we are also including the case where the stray light is not produced 
in the penumbra but further out. Finally, 

C~0.3, (A8) 

a figure coming from the continuum intensity in our umbra when a = 0.05, and when the stray-light intensity is given by the quiet 
Sun intensity. (We use equations IIA2I and IIA4I at continuum wavelengths to evaluate C.) Equations (IA5b (IA6b (IA7b and (IA8b 
yield, 

|A„|<45ms-', (A9) 

which sets an upper limit to the effect of the stray-light on our absolute wavelength scale. A final comment is in order. We do not 
consider the spatial stray-light produced by the spectrograph because it is believed to be negligible. The slit of the spectrograph 
selects only a very small portion of the solar surface. 

INTENSITY- VELOCITY CORRELATION IN LOW RESOLUTION GRANULATION IMAGES 

The observed correlation between intensity and Doppler shift is not one-to-one. As we mention in the main text, a clear 
correlation can be washed out due to the still insufficient resolution of the observations, which may not allow us to see individual 
penumbral convective cells. In order to illustrate the effect, we have carried out an analysis similar to that yielding Figure |5]but 
using low spatial resolution quiet Sun granulation. Figure [TTk shows intensities and velocities obtained from the disk center 1" 
resolution Fe I 15648 A observations described by ISanchez Almeida et al.l (l2003h and lDomfnguez Cerdeiia et"an (l2006h . (These 
spectra were used for convenience, but the behavior should be representative of any other line.) The Doppler shifts have been 
measured as the barycenter of the spectral line, which is fairly symmetric. One can readily see the correlation between bright 
features and upflows, the latter shown as blue contours in the figure. Note how the correlation is not perfect, so that sometimes 
redshift contours overlay locally bright features, and vice-versa. Figure [T7b shows the maps in Figure [T7h smeared with a 
Gaussian 2" wide. The correlation is still visible but it worsens. The 2" resolution maps could be comparable to our penumbral 
data assuming the penumbral convective cells to be 0"1 wide or half our spatial resolution (§|2]l. The granulation cells are some 
1" wide and so half the resolution of FigurefTTb. This simple numerical experiment shows how the quiet Sun convection observed 
with insufficient angular resolution produces maps of intensity and velocity showing only a moderat e correlation, despite the fact 
that the intrinsic local correlation must be very high (see, e.g.. Fig. 3 of the realistic simulations by IStein & Nordlundll998l) . 
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Fig. 17. — Continuum intensity and Doppler velocity contours in a quiet Sun region at the solar disk center. Left: original 1" resolution observations. The 
contours correspond to ±250 m s"'. Right: same observation smeared with a 2" Gaussian, and with velocity contours at ±150 m s"' Note the correlation 
between upflows (blue contours) and bright features. The correlation is not perfect, decreasing as the resolution worsens. 
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